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Abstract 

Q"^ , Some exact solutions for the Einstein field equations corresponding to 

inhomogeneous G2 cosmologies with an exponential-potential scalar 
O^' field which generalize solutions obtained previously are considered. 

I , Several particular cases are studied and the properties related to gener- 

alized inflation and asymptotic behaviour of the models are discussed. 
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I Introduction 



It is generally believed that the scalar field with different kinds of poten- 
tials or different couplings to the curvature makes the space-time inflate and 
isotropize. In this sense inflation is thought as a sort of mechanism that 
achieves the Misner's program [Q, i.e. to explain the large-scale homogene- 
ity of the Universe without imposing very special conditions on the initial 
expansion. This belief is based on the cosmic no-hair theorem proved by 
Wald 0, which asserts that all initially expanding Bianchi models, except 
type IX, with a positive cosmological constant evolve towards the de Sitter 
solution. The same was later proved for a particular class of inhomogeneous 
cosmological models by Jensen and Stein-Schabes 0. It has been pointed 
out @], however, that these theorems exclude a large class of cosmological 
models. In addition, Heusler [^j extended the CoUins and Hawking analysis 
[P] to scalar fleld cosmologies and showed that Bianchi models, except type 
IX, coupled with a scalar fleld with a convex positive potential isotropize only 
if they admit a FRW model as a particular case. Actually the later result 
conflrms that the presence of a scalar fleld does not imply that the model 
isotropizes, and more assumptions about either the form of the potential or 
the coupling between the fleld and the geometry should be then explored. 

Recently, there has been a considerable interest in scalar flelds with the 
Liouville form (exponential) potential. It arises as an effective potential in 
many theories such as Jordan-Brans-Dicke theory, Salam-Sezgin theories, the 
superstring theories and others . On the other hand there exist possibilities 
to flnd exact solutions of the Einstein equations for this potential in a variety 
of cases [H, 0. Moreover, it has been shown that the Einstein fleld equations 
for homogeneous models with scalar fleld with exponential potential can be 
decoupled and, therefore, are reduced to an autonomous system of flrst-order 
ordinary differential equations which can be analysed qualitatively |T^, with 



the de Sitter solution acting as an attractor [11 1, |12|, for particular values 
of the slope of the potential. 

In the inhomogeneous case, little has been done on exact solutions and, 
due to the dependence of the metric on the spatial coordinates, the Ein- 
stein equations can not be reduced to an autonomous system and, therefore, 
the asymptotic behaviour can not be analysed in general. Although, some 
numerical and qualitative results suggest that under certain conditions, in- 
homogeneities can be smoothed and spacetimes isotropize, inhomogeneous 



2 



exact solutions found recently behave in a different way |p. These models 
represent an inhomogeneous generalization of Bianchi I cosmologies, where 
the inhomogeneities are induced either by the irregularities of the scalar field 
or by the initial inhomogeneities in the geometry, due for example to the pres- 
ence of strong gravitational waves. The behaviour of the models suggests that 
more work on inhomogeneous solutions with scalar field and, specifically, on 
exact solutions is needed to clarify the dependence of the inflation and the 
isotropization on the inhomogeneity of the spacetime. 

It has been shown by Wainwright [|l3l that the self-similar solutions serve 
as attractors in the homogeneous case, therefore playing an important role in 
studying the asymptotic behaviour of the cosmological models. The inhomo- 
geneous models found in were self-similar. Thus, it will be interesting to 
see as to whether the models of serve as attractors for the more general in- 
homogeneous solutions, generalizing the property found by Wainwright into 
inhomogeneous cosmologies. 

In this paper we generalize the results of ^ finding a large class of inho- 
mogeneous solutions with scalar field with exponential potential. In Section 
II we describe the geometry and the matter content of the model, and we 
present the corresponding Einstein equations. In Section III we deal with 
these equations and obtain some particular solutions, which can be classified 
into two main groups. We also make a brief discussion about the existence of 
inflation in those cases. We conclude and summarize our results in Section 
IV. 



II The Einstein Equations 

We will consider solutions with one dimensional inhomogeneity, which can 



be described by the generalized Einstein-Rosen space-times admitting 
an Abelian group of isometries G2 and include the Bianchi models of type 
I- VII as particular cases and, therefore, the fiat and open FRW solutions. 
The line element is 

ds^ = ef{-dt^ + dz^) + gabdx''dx\ a,b = 1, 2, (1) 

and the functions / and Qab depend on t and z. 

Assuming for simplicity that the two Killing vectors are hypersurface 
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orthogonal, the hne element (1) may be put in the following diagonal form 

ds^ = ef^'''\-de + dz^) + Git, z){e^^'''^dx^ + e-^^'^'Uy^). (2) 

This line element generalizes one studied in |^ by allowing the transitivity 
surface area G to depend on t and z coordinates. 

The matter source of the metric is that of a minimally coupled scalar field 
with the Liouville type of potential 

V{(j)) = Ae'"''. (3) 

The energy-momentum tensor for the scalar field is 

Tap = (P,a(l),f3 - gal3 ^0,70'^ + '^(0)) • (4) 

This tensor can be rewritten in the form of a perfect fluid stress-energy tensor 

Tafi = (P + p)UaUf3 + pgaf3, (5) 

where Ua, p and p are expressed in terms of the scalar field as follows 

Ua = , , (6) 

P = -^0,a0'" + ^(0), (7) 

P=-^0,a0'"-n0), (8) 

as long as the gradient of the scalar field is timelike. This interpretation will 
be useful to study the kinematical behaviour of the solutions. 

The kinematical quantities relevant to the study of these models are the 
expansion 

e=U^;u9'"' (9) 

and the deceleration parameter 

q = -36-2 (Q^^^a ^ 1q2^ (^q) 



4 



The deceleration parameter will act as an indicator of the existence of infla- 
tion. If g > the solution decelerates in the "standard" way, while q < 
indicates inflation. 

The Einstein equations for the line element (2) and the matter described 
by the stress-energy tensor (4), along with the Klein-Gordon equation for the 
scalar field, are given by 

C C £ r)V 

+ .0. (11) 



G G" 
G~~G 



2efv, (12) 



h-h" + ^h- %h' = 0, (13) 
G G 

where Eq.fpTD is the Klein-Gordon equation, and primes and dots represent 
the derivatives with respect to z and t, respectively. 

We must note that Eq.(12) has the form of an inhomogeneous wave equa- 
tion. On the other hand, the equation for the transversal (or wavelike) part 
of the metric, h, Eq.(13), is very similar to Eq.(ll) for the scalar field. This 
will permit us to solve the Klein-Gordon equation in terms of the metric 
function h. Finally, Eqs.(14) and (15) are the equations for the longitudinal 
part of the metric function /. 



Ill Solving the Equations 

Without any loss of generality, we can write the scalar field in the form: 

= -^log(G) + $ (16) 
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and substitute this expression into Eq.(ll). Then, using Eq.(12) and the 
expression (3) for the potential, the following equation for the function $ 
results: 

$ - $" + ^$ - = 0. (17) 

Gr Gr 

This is the same equation as for the metric function h, Eq.(13), and allows 
us to write the following expression for the scalar field: 

k 

(f) = — log(G) + mh, m = constant. (18) 

It is convenient now to see Eq.(12) not as a differential equation for G, 
but rather as an equation which defines the function /: 

f{t,z) = -k<P + hg(^^-^^ -log(2A). (19) 

We have used by now Eqs.(ll) and (12) and are left with (13)-(15) to 
determine the functions G{t, z) and h{t, z). 

We assume now that the metric function G is separable, as follows: 

G{t,z) =T{t)Z{z). (20) 

Due to the symmetry of the Eqs.(13)-(15) under t-z interchange, we assume 
two different behaviours for function h: 

i) 

e\t,z) = Q{t)Z{zr (21) 

e^{t,z)=T{tyP{z) (22) 

The second ansatz is somewhat similar to that one used by Ruiz and 
Senovilla in their study of inhomogeneous perfect-fiuid G2 cosmologies, [P^], 
while the first one is allowed by the peculiar behaviour of the system of 
differential equations. 

For the first case, Eq.(13) takes this form: 



6 



while for the second case. 



T _ Z'P' P" P' 



nsa 



0,±1), 



(24) 



We wiU study both types of solutions in the following sections. 



III.l Case i) e'^(i, z) = Q{t)Z{zY 

We are now able to write down the form of the desired solutions, in terms of 
T{t), Q{t) and Z{z), which is 

' G{t,z) ^ TZ, 

h{t,z) = \og{Q)+n\og{Z), 

(l>{t,z) = -|log(TZ)+m/i, (25) 

fit,z) = -A;0 + log(f-f)-log(2A), 

V{(f)) = Ae^*^. 

The spatial part of Eq.(23) results in the following three different types 
of solutions: 



Acosh.(az) + Bsmh.{az) if s — 1, 
Z{z) ^\ Az + B if £ = 0, 

Acos{az) -\- Bsn\{az) if e — —1. 



(26) 



On the other hand, if we take the temporal part of Eq.(23), we can 
formally write a relation between the two functions, Q{t) and T{t), as follows: 



InQ(t) 



nea 



T{T)dT 



dt. 



(27) 



This expression, however, is not particularly useful in order to find exact 
solutions, because it yields two intcgro-differential equations. A more prac- 
tical way of handling the equations would be to differentiate the function / 
in Eq.(19) and then substitute / and /' in Eqs.(14) and (15). This will be 
done below. 

We can now easily obtain some particular solutions just by choosing one 
of the three possible cases from (26), and then finding the corresponding 
solution for T{t). 
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III. 1.1 Particular Solutions 



The Linear Case 



The simplest type of solution is obtained when the spatial dependence of 
the metric functions is linear. Rescaling constants A and B we can assume, 
without loss of generality, that Z{z) = z. It can be shown then that the 
function T is given by T[t) = exp (^At + B^, where A and B are arbitrary 
constants. The line element then takes the following form: 



2A 



z - - (-de + dz^) + + z-'^dy' 



(28) 



In the special case when -n?' — 1 the metric represents a homogeneous isotropic 
universe. 

We have looked at some invariant quantities like R^'^Rab-, R"'^'^'^Rabcd and 
have not found any singular behaviour apart from the irregular behaviour at 
z — 0. This singularity disappears in the special case with — 1, when the 
metric is of FRW type. 

Making use of expressions (7) and (8), we can readily calculate the ex- 
pressions for the energy density and the pressure, resulting 

Az-^+V2=:^ f-3 + n2 + 2v/2^ + 3z2) 
^ ^-^^ (2^) 



and 



A ^-3+V2^ (_3 + ^2 ^ 2 V2" 



— z"^ 



P = ^ —t -■ (30) 



The expressions for the expansion and the deceleration parameter are quite 
lengthy, and we will not write them down explicitly here. We will discuss the 
behaviour of the deceleration parameter at the end of this section. 
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The Hyperbolic Case 



The second choice for the function Z is, in the most general form, 

Z{z) = Asmh{az) +Bcosh{az). (31) 

To simphfy the equations, we will take the special case Z(z) — sinh(a2;). The 
general study of the resulting equations for T(t) is more complicated than 
in the previous case, and it is given in the Appendix. There, we will show 
that it is impossible to find an explicit general solution for the system of 
Eqs.(ll)-(15), when we suppose Z{z) — Asmh.{az). In this subsection we 
will present one of the possible particular solutions. 

One can suppose, for instance, that the solution for T could be of the 
same form as for Z, that is, T(t) = Asinh{bt) + Bcosh{bt). If we introduce 
this solution into the equations, we come to the conclusion that A must be 
equal to B, so the actual solution is again an exponential. The complete 
form is 



f Git, z) 

h{t, z) 
<P{t,z) 

fit, z) 

V(0) 



sinh(a2;) exp(6i), 

n log(sinh(a2;)) + 



na 



Ci log(sinh(a^)) + C2t, 

—kCi log(sinh(a2;)) — kC2t + log 

A exp{kC2t) [sinh(a2;)]'^'"^ , 



(32) 



2A 



(fc^-2n^+2)(V^)-fc 
2 



with Ci 

1 ^ _ a^(fc^-2n2+2)(V^)-fcb^ 

The expressions for the energy density and pressure, Eqs.(7) and (8), 
become more complicated for this solution: 



Ke-ycp{kC2t) sinh(^^i)(a^) ( '^1 cot'(a^) + 1 



62 -a2 



p = Aexp(A;C2t) sinh^^^^^) [az] 



' CI- cy 

62 



cot^(a2;) - 1 



(33) 
(34) 



The corresponding line element is: 



fo2 



-kC2t , 



sinh-''^' {az){-df + dz') + 



2A 



+ ( sinh^+"(a^)e 



dx^ + sinh^ ''{az)e 



(35) 



From the expressions (33) and (34) we deduce that there is a periodic 
singular behaviour in z. 



The Trigonometric Case 

This is the last type of solutions, which is not very different from the 
previous one. We can try to find again a particular solution with an expo- 
nential behaviour in time, as it seems to be the most common among all our 
solutions. In fact, there is a particular solution of this kind, namely: 



G{t,z) 

h{t,z) 
</){t,z) 

f{t,z) 

V{<t>) 



— sin(a2;) exp(6t), 



na 



n log(sin(a2;)) — t, 



Ki log(sin(a2;)) — K2t, 

—kKi log(sin(a2;)) + kK2t + log 

A exp{—kK2t) [sin(a2;)]'^''^^ , 



2A 



(36) 



/ 



1 (A;2_2n2+2)(l/2)_fc 

where Ki — 



V 



a2(fc^-2n2+2)(i/2)+fcb2 



26 



As can be easily seen, this solution is very similar to the hyperbolic one. 



III.1.2 Do The Models Inflate? 

As to the existence of inflation in these models, we have to say that all the 
expressions for the corresponding deceleration parameters are quite compli- 
cated, and therefore, little information can be obtained from their explicit 
form. However, we can look at the behaviour of q against one of the variables 
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t or z for some particular values of constants. It can be shown analytically 
that, in all of the cases discussed above, the sign of the deceleration param- 
eter does not depend on time. This means that the inflation has a spatial 
character. Some of the regions, depending on their location in the universe, 
will inflate forever, but some will never experience the inflation. 

The most remarkable feature is that all of the models wc have looked at 
show a strong dependence on the spatial coordinate z. Moreover, given the 
same kind of solution with shghtly different values for the constants, we get 
a quite different behaviour of g as a function of z. In the case of periodic 
solutions {Z — sin(a2;)) the inflating and non-inflating regions interchange 
periodically. 

When a cosmological model deviates strongly from homogeneity it can 
be extremely difficult to see whether the model inflates or not. When this 
happens, however, we can use a somewhat weaker way to deduce about the 
inflationary behaviour of the model, just by looking at the fulfillment of 
the strong energy condition. The breaking of this condition is known to be a 
necessary for inflation to take place. Therefore, we can study the behaviour of 
the quantity £ — p+3p for different values of the variables and the constants. 
When this quantity is such that S <0 the strong energy condition is violated, 
and inflation might take place. 

If we do so for the first of our models, the linear case, using expressions 
(29) and (30), we can see that S is always negative (again, the sign of this 
magnitude has no dependence on time, since this coordinate always appears 
in the expressions of p and p in the exponential form). Therefore, the model 
can infiate. For the special case with = 1 we have £ — 0, which corre- 
sponds to an equation of state of the form p = — p/3. 

In the hyperbolic case, making use of Eqs.(33) and (34), we can see that 
the strong energy condition is broken only for some values of the coordinate 
z and depend strongly on k. The energy condition is always broken for small 
z. \{ ^ 2 the energy condition is broken for a wide range of values of 
z. When k grows, this range decreases strongly. Therefore, we can say that 
inflation is always possible near the hypersurface 2; = 0, but is less probable 
in other regions. 
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III.2 Case ii) 



e''{t,z) =T{tYZ{z) 



As we have discussed at the beginning of this section, we can use a different 
assumption for the form of the metric function /i, using Eq.(22) instead of 
Eq.(21). If we do so, we must deal with Eq.(24) rather than Eq.(23): 



T Z'P' P" P' 

^f~~zp'^y~^ 



nea 



0,±1), 



and again we can easily solve the temporal part of this equation, which is 
equivalent to Eq.(26), but in this case with t as the variable: 



Acosh.{at) + Bsmh{at) if e — 1, 
At + B if £ = 0, 

Acos{at) + Bsm{at) if e — —1. 



(37) 



Once more we have three different cases. 

Similarly to the previous case, we can write down the general form of the 
solutions, by taking G from Eq.(20), h from Eq.(22), and (f), f and V from 
Eqs.(18), (19) and (3) respectively. The most general solution can be written 
then in the following form: 



G{t,z) 

h{t,z) 

V{<l>) 



TZ, 

nlog(T)+log(P), 
-|log(rZ)+m/i, 

-fc + log — — 



(38) 



log(2A), 



Again, we can try to find the most general solution corresponding to each 
particular case of T{t), or at least, find some particular solutions that could 
be of interest in each case. 
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III. 2.1 Particular solutions 



The Linear Case 



If we take the simplest case from Eq.(37), which is the hnear solution for 
T{t), we obtain a general solution analogous to Eq.(28), that is, 

ds^ = te^{-df + dz^ + t^dx^ + r^dy^). (39) 
The expressions for the energy density and the pressure are now 

p = ^e-^ {t-' - 3t-') , (40) 

and 

p= V(r^ + r^). (41) 



This is an interesting solution owing to its simplicity, which allows us to 
calculate simple expressions for the expansion and the deceleration parame- 
ter: 

2eft6(i-3_i-i)2 

,^e-- (3t^-i)(t^ + 9) 

Looking at the sign of g, we can now easily see that inflation starts at 
This universe expands in the standard way before t = l/\/3 and then inflates. 



13 



The Hyperbolic Case 



If we consider the second case, T(t) oc sinh(at), we can take a particular 
solution which is the equivalent of equation (32), interchanging variables. 
The spatial part of function G{t, z) is an exponential, as in the previous case. 
All the particular solutions found appear to have an exponential behaviour 
in z. The solution is the following: 



f G{t,z) 
h{t,z) 

<P{t,z) 
f(t,z) 

y(0) 



sinh(at) exp{bz), 
log(sinh(ai)) + —z, 

— A;log(sinh(at)) — — 

A;^log(sinh(ai)) + — 

Ae^*^, 



z + log 



(45) 



2A 



with b = a\ 



The expressions for the energy density and the pressure are as follows: 



P 



P 



(sinh(ai))'=^ 



(sinh(ai))^^ 



4A 

a^-b^ 



k^a^coth^at) - — 



4A 



k^a^coth\at) - — 



e fa^ + b' 



+ A 



A 



(46) 



(47) 



It is worthwhile to mention that among these solutions we have also found 
a solution which represents a universe with a cosmological constant and a 
vanishing scalar field. The corresponding line-element reads as follows 



ds'^ = J (-dt^ + dz^^ + smh^{t)dx'^ + cos^ {z)dy'^ 
This line-element can be put into the following form: 



ds' 



1 
A 



-dt^ ^ Ks\W{t)dx'^^ + ^ (dz^ ^ Kcos^ {z)dy 



A 



(48) 



') 
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which can be recognized as a Kantowski-Sachs model. These models have a 
G4 symmetry which acts multiply-transitively. Although this solution does 
not appear in the literature we assume that it must have been derived by 



Kantowski in Ref. [ICl, as mentioned in Ref. 17 



In all these solutions the sign of the deceleration parameter is spatially 
homogeneous, this is to say that it does not depend on the spatial coordinate 
z. The results of this subsection can be easily generalized to the case of 
trigonometric functions oc sin(at)). In this case one can speculate about 
the possibility of the periodic in time inflationary behaviour - the universe 
in which the inflationary and non-inflationary periods alternate each other 
in time. 



IV Qualitative Asymptotical Behaviour of the 
General Solution 

We will try to solve the general system of equations given by (11)-(15) with 
the restrictions (18), (19), (20) and (21). As we have seen in Section III, 
this leads to Eq.(23), from which the spatial part of the solution is Eq.(26). 
The general solution is then given by Eq.(25). We have now three equations, 
which are Eqs.(13)-(15). Eq.(13) becomes Eq.(23), as explained above, and 
the new expressions for Eqs.(14) and (15), after introducing the general form 
of the solution (25), follows: 



^ Z' f 



T 



[1 + k'^ — kmn) — f + {2m'^n — mk + n 



Q 
Q 



(49) 



Z^ 

~ T 

T 

+ 



1 



-n 



1 

— V m n H kmn 

2 4 



+ / 



— km 



TQ 



Z^ 

'e_l\T^_.T fl ,\Q^ 

T ■ V4 2;t2 V2 y g2 j^q- 

We note that these equations are not completely separable, since the deriva- 
tives of the function / depend both on t and z, which can be seen from the 
following expression: 



(50) 



fit,z) 



k' 



log(TZ) — kmn\og{Z) — km\og{Q) + log 



z 



2A 



(51) 
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It can be shown, however, that for the three special cases under consideration, 
/ and /' are functions of t and z alone, respectively. 

Separating Eqs.(49) and (50) into a spatial and a time-dependent part 
and introducing the function / into the time-dependent part of the Eq.(50), 
we finally get 



— + -^+ - + mM^-( 2m^n + n) --^ = a, 52 



where a is a separation constant. Introducing now two new variables 

« = ^ and ^ = ^' 
the pertinent equations (23) and (52) can be written as: 



(3' + a(3 = nea?, 

a' - Ka^ + - 2Mnap - a = 0, 



(54) 



where 

Pi 1 

K = — + - and M = - + ml (55) 

4 2 2 ^ ^ 

Expressing a from the first equation and substituting it into the second, 

we obtain a single equation for (3. Introducing then the new variable y{l3) = 

f3', we get a non-linear differential equation of the form 

y'y = W) + + W)y\ (56) 

where are certain non- linear functions of (3. This is a so-called Abel 

equation of the second type [0, whose solution can not be inverted to obtain 



/3 as a function of y. Therefore, it is not possible to get the desired general 
solution for the system of Eqs.(49), (50). However, we can try to make an 
analytical study of the system (54) and get some indications of how the 
behaviour of the solutions is. 

To achieve this, we have studied the phase space of the system (54), which 
is equivalent to that of Eq.(56). The system has two equilibrium points which 
are saddle points, therefore unstable. These two fixed points correspond to 
inhomogeneous solutions satisfying Q(t) = e^"* and T{t) = e"""*, with ao and 
Po being the coordinates of these points in the phase space. In Figures 1 
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and 2 we represent the phase space of the system (54) for different values 
of the constants. The trajectories of the solutions have a slightly different 
behaviour in both cases, since in Figure 2 the two saddle points are closer to 
the /? = axis. In both figures, we have assumed that e = 1 in (54). In the 
case that e = and cr = 0, which corresponds to solutions linear in z, there 
is a special family of stable solutions, those lying on the negative part of the 
a-axis, that is, with /3o = and ao < 0. All these solutions evolve towards 
a fixed point at the origin, of a static nature. This behaviour is shown in 
Figure 3. 

The same study can be done interchanging t and z coordinates, obtaining 
the same system (54) as the final result, but this time a = ^ and P = 

and we have a system of equations depending exclusively on the spatial 
coordinate. 

V Conclusions 

We have studied a large class of G2 inhomogeneous universes with a scalar 
field with a potential of the Liouville type. After writing the Einstein equa- 
tions we have seen that there are two main classes of solutions. In the first 
case, we have supposed that all the metric functions have a similar spatial 
dependence. As a result of this assumption, we obtain that the function 
Z(z), which gives this spatial dependence, must have the form of one of the 
three cases stated in Eq.(26). Another consequence is that the sign of the 
deceleration parameter, which informs us about the existence of inflation in 
these models depends only on spatial coordinate z, that is, some parts of 
the universe will inflate while the other will not, at any time. On the other 
hand, the second type of solutions represents somewhat an opposite case: all 
the metric functions present a similar dependence in time, and the function 
T(t), which drives this dependence, is given by Eq.(37). The sign of the 
deceleration parameter depends now on time only, so that the whole of the 
universe inflate at the same time, but the periods of the inflation and the 
duration depend on the initial conditions. 

We have looked at the asymptotic behaviour of the solutions at late times 
studying the system of equations (54). We have found no stable equilibrium 
points for this system. 

This work mainly was concerned with exact solutions. We have seen that 
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by introducing the large degree of inhomogeneity the cosmological model can 
have a very complicated structure. 
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FIGURE CAPTIONS 



FIGURE 1: Phase space corresponding to the system of equations (54) and 
the following values for the constants: n — 1, a — 1, k — ^/2, m — V2/2 and 
a — 1. e is equal to 1, so we are in the hyperbohc case of section III. 1.1. 

FIGURE 2: Phase space corresponding to the same system (54), but for 
these values of the constants: n — 1/2, a — 1, k — 2, m — 1 and a — 0. e is 
still 1. 



FIGURE 3: Phase space corresponding to £ 0, that is, to the linear case of 
section III. 1.1. There is a collision of saddle points at the origin. The values 
for the constants are those of Figure 1, but with a — 0. 



